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We study the dynamics of an initially inverted atom in a semi-infinite waveguide, in the presence 
of a single propagating photon. We show that atomic relaxation is enhanced by a factor 2, leading 
to maximal bunching in the output field. This optimal irreversible stimulated emission is a novel 
phenomenon that can be observed with state of the art solid state atoms and waveguides. When 
the atom interacts with two one-dimensional electromagnetic environments, preferential emission in 
the stimulated field can be exploited to efficiently amplify a classical or a quantum state. 
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I. INTRODUCTION 

Stimulated emission is a central concept of quantum optics. Initially introduced by Einstein [1 in the case of a 
macroscopic number of emitters and electromagnetic modes, it describes the increase of light-matter coupling with the 
number of photons in the mode. Because of this non-linear behavior, an initially inverted gain medium preferentially 
emits light in the stimulated mode, a feature that led to the development of lasers and amplifiers technology [2J. 
Latter experimental developments have allowed to design optical media sensitive at the single photon level. One 
paradigmatic example is given by a two- level atom strongly coupled to a monomode cavity [3j|4]. In this specific case, 
spontaneous emission is reversible and gives rise to the so-called vacuum Rabi oscillation, whereas the presence of a 
single stimulating photon enhances the frequency of the oscillation by a factor y/2 [5] . More recently, another class of 
giant optical non-linear medium has emerged, like atoms weakly coupled to directional cavities [6j [7] , superconducting 
qubits in circuit QED [8 , or quantum dots in photonic wires [9]. Such atoms that interact with only one direction 
of the light field open promising perspectives in quantum communication and information processing, as they are 
expected to provide efficient single photon transistors [lOj [11] or photonic gates [I2j [13]. These one-dimensional 
systems evidence new fundamental effects which have been extensively studied both experimentally p~4[ [15] and 
theoretically [T6UT8] , just to mention a few examples. 

Here we explore the dynamics of stimulated emission for a two- level atom embedded in a semi-infinite waveguide. 
We show that a single resonant photon, properly shaped, propagating in the waveguide can shorten the atomic 
lifetime by a factor 2, leading to significant bunching in the output light field. Such optimal irreversible stimulated 
emission is a new phenomenon: in the context of cavity quantum electrodynamics, stimulation by a single-photon is 
either optimal, but reversible or irreversible, but not optimal. Building on this effect, we revisit the case of an atom 
interacting with two one-dimensional electromagnetic fields, as it is the case for a two-level atom in a transmitting 
waveguide, or a lambda shaped three level atom in a half waveguide. Under optimal conditions, light emission is twice 
more probable in the stimulated field than in the field prepared in the vacuum state. This effect is discussed in the 
context of amplification and quantum cloning. 

II. MODEL AND METHODS 

The system under study is represented in Fig{l] An initially inverted atom is embedded in a waveguide where only 
one direction of propagation is allowed, which could be realized by a semi-infinite waveguide, closed by a perfect mirror. 
Atomic emission may be stimulated by a single photon pulse propagating in the waveguide. The total Hamiltonian is 
given by 

H = i^iD + H atom + i^ int , (1) 
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which is the same system analyzed in Ref. [19 j. It describes an emitter interacting with a one-dimensional electro- 
magnetic field, given by the Hamiltonian Hm = Yl^Lo ^ v a t a vi where is is the frequency of the a v mode [T9] . The 
atom with frequency is a is modeled by H atom = his a cr+a_, where cr_ = \g)(e\ and <j + = al. The dipole interaction 
is given in the rotating-wave approximation by 

oo 

Him = -i ^2 hg » l a v a + ~ Hx -] ( 2 ) 
for the atom centered at the origin, where g v is the coupling strength. 
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Figure 1: Scheme of the two- level atom stimulated by a single photon pulse of linewidth A. The spontaneous decay rate is Y 
and the transition frequency is is a- The direction indicated by the arrows is the only allowed one, which models the presence 
of a mirror close to the atom, at position ta- The light field is detected by a photodetector positioned at rd, put sufficiently 
far from the atom. 



The solution of the problem is Hamiltonian, so the total number of excitations is conserved. An Ansatz is constructed 
for the complete wavefunction \ip(t)), 

CO CO 

|V(*)> = 4|e,0)+ ^."'(*) 44'M>, (3) 

that evolves according to the Schroedinger equation, 



CO 



d = -i("A + y)^v - 2 ^2 dvfoy, (4) 



and 



dt 

v'=0 



A. Excited-state amplitude in real-space representation 

For the following analysis we adopt a real space representation of the quantum state, ip(r,t) = Y^Lo^i^) e lkl " r , 
where k u = \k u \ = is / c. With the help of the transformation Y v v ip u e lkljr = —icd r ifj u e tk " r , we rewrite Eq. Q in 
the form [J^ + c ~§p] ^( r ? ^) = ~^ v a ^( r ? — 2 v , g v > § v ye %kur '. In the trivial uncoupled case g v = 0, for instance, the 
solution to the equation above is simply a product of a free time-dependent atomic evolution f(t) = ex.p(—iis A t) and a 
propagating pulse p(r — cb), where p(r) is the initial pulse wavefunction. In this paper we integrate Eq.([5| choosing the 
initial condition (j) v y{G) = 0, so that the atom is initially in the excited state with a single-photon propagating in its 
direction. As usual, we assume the fast rotating reference frame, ip(r : t) = ifj(r, t) exp(—iis A t) ex.-p(—iis L (t — r/c)), where 
is l is the central frequency of the incident wavepacket. Two terms in our equations deserve further considerations. The 
first is Jq dt' 'ijj(r — c(t — t') , t') x gl,e~^ u '- v A){t-t') wn i cn equals to (T/2) ^(r, i), under Markovian approximation. 
The decay constant is T = 27rgl,5(is f — is A ), as obtained by the Fermi golden rule for spontaneous emission. The 
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other term can be shown to satisfy TQ(r)Q(t — r / ' c)e~ %bLl "> /c ijj(—r,t — r/c), under reasonable approximations, namely, 
9v ~ 9v A and the continuum limit J2 U — >> j dvp\D- The detuning is defined as Sl = — and the ID density of 
modes is puj = L/(ttc). Finally, we have been able to eliminate self-consistently the dependence on the two-photon 
amplitude, finding 



d_ d_ 
dt dr 



rj>(r,t) = ~ ip(r,t) - T 0(r)G(t - r/c)e- i5Lr/c tP(-r,t - r/c). (6) 



Note that in the high detuning limit, Sl ^> T, the solution is simply ijj(r,t) = e ^ t ijj(r — ct, 0), i.e., the product of 
a decaying atom and a freely propagating photon. The general solution for Eq.(|6| reads 



$(r,t) = ^(r-ct,0)e-^-r6(r)e(t-r/c)e-^e-(S-^) (£ - r/c) / e (*- iSL ) t ' $(-ct f , 0) dt\ (7) 

Jt-r/c 

where the initial condition of the wavepacket is written in ^(r, 0). In the following we restrict the study to the 
case of an incident photon of exponential shape, as if it had been spontaneously emitted by a neighboring atom of 

natural linewidth A and central frequency vl-> i.e., ^(r, 0) = A/"e( T+zl/£y )^6(— r), where M 2 — 27rpi£>A stands for 
normalization. To obtain inversion of population of a two-level system in an experiment, a pulsed-laser excitation 
can be used [9]. Another option is to incoherently pump the two- level atom through a third level, which is properly 
described by a different formalism [20 . 

B. Two-photon amplitude in real-space representation 

The two-photon component are also defined in real space as </>( r i? r 2^) = v ' 4>vy(t)e zk " ri+zku ' r2 ' Within the 
same approximations done before, we find 



0(ri, r 2 , t) = sj [6(t - r 2 /c)6(r 2 )^(ri - r 2 , t - r 2 /c) + Q(t - n/c)e(ri)^(r 2 -r u t- n/c)] (8) 

written in the original reference frame. In the limit of vanishing interaction between the incoming photon and the atom 
(e.g., Sl ^ T), the two-photon wavefunction can be shown to consist in a symmetrized product of two independent 
single-photon wavefunctions, one describing spontaneous emission (^sp.em.fV^) = exp[— (T/2 + ivA)(t — r/c)]) and 
another describing the free propagation of the field (^free prop.( r — c£, 0)). In the limit of ^> the excitations 
are entirely in the field and the dynamics also reduces to free propagation, so that one can define the function (j)^ 
checking 

0(^1,^2, too) = 0oo(ctoo - ri,ctoo - r 2 ). (9) 

III. TIME RESOLVED SIGNATURES OF STIMULATED EMISSION 

A. Excited state population dynamics 

We first investigate the dynamical signature of stimulated emission and compute the excited state population 
Pee(t) = (e|Trfi e id[|V ; (t))(V ; (t)|] \ e ) • Given our choice for the initial state, 



Pee(t) 



" (r+A)t jl + |1 + Q\\e At - 1) + |Q| 2 ^(e F£ - 1) - 2Re 



,(10) 



where we have defined the factor Q = 2T/(r — A — 2iSl) and Re stands for the real part. First, we note that in 
the strongly detuned case Sl ^> I\ the atom is transparent to the incident photon, so that p e e(t)\s L -^oo — exp(— IY), 
which is spontaneous emission. We then focus on the resonant case and study the influence of the width of the pulse 
on the dynamics of the atomic relaxation. A convenient signature is given by the adimensional effective lifetime 
T e ff = T J °° p ee (t) dt, which reads 
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and is plotted in Fig. [2j The case A ^> T corresponds to a very short pulse in time and also gives rise to free 
propagation of the pulse followed by spontaneous atomic decay (r e ff = 1) as the spectral overlap with the atom is 
negligible. On the contrary, a highly monochromatic photon (A <C T) increases the effective lifetime due to the 
efficient absorption of the incident wavepacket, after the atom has spontaneously relaxed. Shortening of the atomic 
lifetime induced by stimulation can be observed for T < A < 100r. If A = T, the integrated effect of the stimulation 
exactly compensates for the total absorption and r e ff = 1. The corresponding dynamics is plotted in the inset of Figj2j 
where a fast decay is followed by a re-excitation induced by the absorption of the tail of the photonic wave packet. 
The case where A = 3r gives the shortest effective atomic lifetime, namely half the spontaneous emission lifetime. In 
this situation, a fast population decay induced by stimulated emission takes place, minimizing the reabsorption effect 
and leading to optimal irreversible stimulated emission. Note that for times t < T _1 (resp. t > T -1 ) the population 
is bigger (resp. smaller) than the reference exp (— 2Tt) so that the effective lifetimes of both curves are equal. 

At this point we mention that an intuitive derivation of the optimal effective lifetime can be obtained by modelling 
the ID atom as an ultimate gain medium, e.g. a single emitter initially inverted, irreversibly decaying into a collection 
of modes {rij}. Einstein rate equations for the excited-state population can be written [TJ|2] 

d ? 

— p ee = -r(l - /3) (1 + Til) pee ~ T(l + n k )j3 p ee , 

where n k is the number of photons in the stimulating mode, n\ the number of photons in the other modes. j3 is the 
fraction of coupling with the one-dimensional (ID) channel with respect to the 3D continuum of modes [20] . Let us 
assume that n\ = 0. In the conventional 3D scenario, f3 <C 1, the emitter's decay is not modified unless > In 
the ID case under consideration, j3 = 1, a single photon (n^ = 1) is enough to stimulate the transition and shorten 
the lifetime by a factor of 2 (p ee = — 2rp ee ). 

Even though it is intuitive, this simplified picture does not capture all the physics of the problem. The atomic 
evolution presented here does not obey a rate equation, justifying the Hamiltonian resolution adopted here. To give 
an example, let us consider the best stimulation condition, A = 3I\ In that case, 

p ee = -2e- 4rt + 3e- 3rt , (12) 



so 

^ = 8 r e - 4rt - 9r e - 3rt + -2r Pee , (13) 

clearly evidencing a difference between standard laser systems and the present scenario. 

Emitters weakly coupled to directional leaky cavities are often suggested as potential one-dimensional atoms O [7] . It 
is worth noticing that atomic emission cannot be stimulated by an additional photon initially prepared in a dissipative 
monomode cavity. As a matter of fact, this photon would escape the cavity in a typical time way too fast to 
stimulate any atomic emission that would take place on n/g 2 , where the atom-cavity coupling strength g checks in 
the weak coupling regime g <C n. Hence, stimulated emission by a single photon can only be simultaneously optimal 
and irreversible with pulse-shaped photons propagating in broadband waveguides. 



B. Temporal correlations in the output field 



A consequence of the atomic relaxation enhancement is the emergence of bunching in the output field. This kind of 
photonic temporal correlation finds application, for instance, in quantum lithography [2TJ [22] . Photodetection signals 
are registered with a detector positioned at ^> c/Y as pictured in Fig [I] so that all the excitations are in the light 
field. This regime corresponds to free propagation, so that the characteristics of the field only depend on the variable 
r — ct and eq.([9| is valid. We shall use the reference frame of the photodetector, of origin and td = r^/c. The 
density of probability to detect one click at time t and one click at time t + r is obtained from the second-order 
correlation function [23] that in our case checks G^ 2 \t^t + r) = |<^oo(c£, c(t + r))| 2 . From Eqs.([7|) and (JS]), we find for 
r > 

G (2) (M + r) = Are-( p+A ^ (l + Q)e-5 r + (l-Q)e-(^+^) r (14) 

where Q = 2T/(T — A — 2i5l) is the effective coupling factor as defined in the atomic population dynamics. A clear 
interpretation can be given to the expression above. Between time t and time t + r, the system is projected on the 
single excitation subspace, giving rise to two possible situations. In one case, the first click comes from the incident 
field. The remaining excitation is in the atom, that will finally spontaneously relax: this corresponds to the dynamics 
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Figure 2: Effective population lifetime r e ff as a function of A. The dashed line indicates the optimal value r op t = 1/2. In the 
inset: excited-state population p ee decay in time. The blue dot-dashed curve, for A = T, shows stimulation for very short times 
and reabsorption after 2T -1 . The red solid curve represents the maximally stimulated relaxation, for A = 3r, which for times 
smaller (bigger) than T _1 goes above (below) the reference exp (— 2Ft) (green dotted). The black dashed curve represents the 
spontaneous emission evolution exp (—Ft). 



exp [— Tr/2] of weight 1 + Q. In the other case, the photon emitted by the atom clicks before the incident one. This 
second situation, that corresponds to stimulated emission, gives rise to the component exp [— (A/2 + iSl)t] weighted 
by 1 — Q. The condition 1 + Q = allows to cancel spontaneous emission and to maximize stimulated emission. This 
is obtained for A = 3T and 6l = 0, namely the very same condition that minimizes the atomic lifetime. 

This optimal regime for stimulated emission leads to the emission of both photons in a typical time 1/3T, giving 
rise to bunching in the output field. The effect can be observed on Fig(3j where we have plotted the probability of 
detecting the two photons within a time Tfi na i (jQ Tfinal dr J °° dtG^ 2 \t^t + r)) as a function of Tfi na i, at resonance, for 
different values of the parameter A. The fastest convergence is obtained for A = 3T. Keeping this optimal value of 
A we have plotted the same function, on Fig. |1J for different detunings, clearly showing as well the importance of the 
resonance on the stimulation efficiency. 




Figure 3: Probability of detecting 2 photons within a time interval Tfi na i as a function of Tfi na i for A = O.ir (blue), IF (green), 
3T (red), and lOOr(black). 
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Figure 4: Probability of detecting 2 photons within a time interval Tfi na i as a function of Tfi na i for an incident photon of different 
detunings: 5l = (red), 5l = 5r (blue), and Sl oo (green). 



IV. POTENTIAL OF ID ATOMS FOR CLASSICAL AND QUANTUM AMPLIFICATION 



Interesting applications of stimulated emission rely on the preferential emission of light in the stimulated mode 
rather than in other empty modes. With this aim, we revisit in this section the case of an atom interacting with two 
one-dimensional electromagnetic fields a and 6, one containing a single propagating photon, and the second one in the 
vacuum state. Two paradigmatic systems are considered in the light of the time-resolved study performed above (see 
FigJEJ: a two- level atom in a transmitting/reflecting infinite waveguide [24], and a lambda shaped three level atom in 
a semi-infinite waveguide closed by a mirror [25]. The fields a and b correspond respectively to photons propagating 
to the right and to the left, or photons of two orthogonal polarizations in the half waveguide. 




Figure 5: Comparison between different types of two-continua of modes in ID atoms. On the left, the reflection ^/transmission 
a are the two degrees of freedom. On the right, the degrees of freedom are represented by the a/b polarizations. The decay of 
the excited-state population is the same in both cases. In contrast, the probabilities of emission in each channel, p a b and pbb, 
crucially differ. 



7 



A. Models and atomic population dynamics 

1. A two-level atom in a transmitting waveguide 
The Hamiltonian for such system reads H = Hm + H atom + Hi nt and the new interaction is given by 

oo 

Hint = -i ^2 hg » + M a + ~ Hx -] ' ( 15 ) 

where a v is the forward and b v is the backward propagating modes [26]. Following the formalism in Ref.[24], one can 
decompose the problem in two parts by using the even e v — (a v + b v )/y/2 and odd d v — (a v — b u )/y/2 field modes. 
The new interaction Hamiltonian depends only on the even modes, 

CO 

H int = [<W - H.c.] , (16) 

while the odd modes contribute only to the free evolution. The spontaneous emission rate V of the transmitting ID 
atom is doubled with respect to the case of an atom in a closed waveguide, i.e. V = 21", which is due to the doubling 
of the resonant modes available for atomic relaxation. 

In the case of stimulated emission under study, the initial state of the field writes |l a ,0&), which checks |l a ,0k) = 
(|l o ,0 e ) + |0 o , l e ))/V% where (resp. |l e )) is the corresponding propagating photon in the odd (resp. even) mode. 
Hence, the system can follow two equiprobable paths, one giving rise to stimulated (resp. spontaneous) emission, so 
that the excited state population of the full transmitting/reflecting ID atom checks 

in terms of the population p|| mi (t) of the two-level ID atom closed by one mirror. 

As previously, the minimal value of the effective atomic lifetime for this system is obtained at A = 3r / '. Because 
of the intrinsic presence of spontaneous emission in the odd modes, the effective lifetime is averaged and increases to 
T Q pt = 3/4 as demonstrated in [24] . 



2. A lambda shaped ID atom 

In the case of a three level atom in the lambda configuration coupled to a half waveguide, the initially excited atom 
\e) can decay by emitting an a or a b polarized photon, ending up in state \g a ) or \g^) respectively. The interaction 
Hamiltonian writes now 



H int = -i^2hg u {a v (j\ + b v a\) - H.c. 



(17) 



where cr a = \g a )(e\ and = |#&)(e| are the lowering atomic operators. As above, the spontaneous emission rate is 
doubled, and the system evolution in the presence of a single propagating photon in mode a also splits into two paths, 
namely spontaneous emission in mode 6, or stimulated emission in mode a. This formal analogy leads to the same 
conditions of minimization of the effective atomic lifetime. The crucial difference between the two setups (transmitting 
waveguide and lambda atom) appears in the characteristics of the fields, as clarified in the next section. 



B. Optimal emission in the stimulated mode — application to amplification 



(2) 

To study the emission in the stimulated mode, a convenient quantity is the density of probability G a a(t,t + r) of 
detecting a click at time t and t + r in that mode, that is given in both systems by 



G ( 2 a ) (t, t + r) = Ar'e-( r '+ A ) t (1 + Q f ) e"^ T + (1-Q f ) e"( * +Ul 



(18) 
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where Qf = V /(V — A — 2z5l) = Q/2 is the new effective coupling factor. Again, the term evolving like e r r / 2 
corresponds to the case where the first click comes from the field, and is followed by a spontaneous emission of the 
atom in mode a. The second term e _ ( ~2+ i5l ) t i s due to stimulated emission of the atom in mode a, followed by a 
second click coming from the field. Extinction of spontaneous emission is realized when 1 + Qf = 0. This yields 
A opt — 2r', a condition which does not minimize the effective lifetime, as it was the case in the half waveguide case. 

A consequence of the optimal correlation is found in the probabilities of photon emission. By suppressing the 
spontaneous emission path, one avoids the atomic emission in the mode b once a first photon a has been detected, 
strongly inducing the emission of light in mode a. Naturally, this condition also maximizes the probability to find 
two photons in the mode a in the end of the interaction, which reaches p aa = 2/3 [24] [25] . This is the maximal 
probability to clone the state a. In the case of the lambda atom, this property can be used to clone the quantum 
state of the incoming photon [25]. A convenient figure of merit for this device is the fidelity T = p aa + p a b/2, which, 
in the optimal point, exactly reaches 5/6, namely, the optimal bound for quantum cloning. Because of the rotational 
invariance of the global system, the cloning is universal, a highly desirable property for state estimation and quantum 
cryptography. 

On the other hand, the transmitting atom can be used as an ultimate gain medium to efficiently amplify the classical 
state of the photon, encoded in the direction of propagation. The convenient criterium in that case is the amplification 
ratio (called visibility in [10] and gain in [11]) A = [T^ pt — T J \Sl — >• oo)]/J rT (^L — >• oo), where the parameter T T 
now represents the transmission fidelity of the system. Remembering that Tmax = 5/6, and J- t (Sl — >• oo) = 3/4, one 
simply finds the maximal amplification ratio A m ax = 1/9. Contrary to the case of the lambda atom, this upper bound 
is not reached here. The transmission fidelity is plotted in fig|6] as a function of the packet linewidth, its maximal 
value being T^ pt — 97,5% x (5/6). This difference comes from the fact that after having spontaneously emitted a 
photon in mode 6, the atom in state g can still interact with the incoming field in mode a. Thus, it is possible to 
finally get two photons in mode 6, so that p^ ^ 0, lowering the fidelity. This situation is forbidden for a lambda 
atom that relaxes in state g^ and becomes transparent to the incoming field. We find A = 1/12 « 8,3%, that means 
almost one order of magnitude higher than previously reported gains, around ~ 1% p~Qj[TT], working in the regime 
where the probe is continuous [20] . 




1000 
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Figure 6: Transmission fidelity T T as a function of the packet width A. The optimal value is T^ vt — 97, 5% x (5/6), found for 
A = 3r x . The dashed line indicates the maximal value T^ax = 5/6. 



In a realistic scenario, two noise sources must be taken into account, namely, the decay rate into the environmental 
3D channels 7 and the pure dephasing rate 7* present in solid-state systems. The former is usually quantified by 
the parameter (3 = T/(T + 7) which can reach 0.98 in ID nanophotonic systems made of photonic wires [9 or ID 
waveguides in photonic crytals [27] [28], and almost 1 in circuit QED [12]. Pure dephasing rates of 7* « O.ir have 
been measured in quantum dots [29] and superconducting qubits [8 . The stimulation depends on the coherence of 
the atom-field interaction, which can be shown [20 to depend on f3 and 7* roughly as ~ — 7*/T), at /3 « 1 and 
7* <C T. One can thus define a factor of trust fa ~ — 7*/r) that estimates how close to the optimal values the 
realistic fidelity can be. For the above mentioned systems, it can be expected of the order of fr ~ 90%. In addition 
to building cleaner systems, dynamical decoupling approaches [30] can be used to reduce dephasing. 
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V. CONCLUSIONS AND PERSPECTIVES 



We have shown the influence of the incoming photon on the atom decay as a function of the packet shape. An 
irreversible and maximally accelerated stimulated emission occurs for the broadband mode-matching condition where 
the incoming photon is three times shorter than the spontaneously emitted one. We have also studied the influence of 
stimulation on the two-photon correlation function, which shows optimal photon bunching. Finally, we added a second 
one-dimensional field to explore quantum and classical amplification. In the quantum case, the possibility to achieve 
universal optimal cloning of polarization has been presented. In the classical case, amplification has been shown in the 
average transmitted field, which reaches 97.5% of the ideal case. This effect has led to a transistor-like amplification 
that can overcome the continuous- wave approach by a factor 8. The photonic propagation makes the reported effects 
specially attractive as far as realistic implementations of quantum information processing are concerned. 
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